{% raw %}
All result for the 8a.1 benchmark specification.
from IPython.display import HTML
HTML('''<script>
code_show=true;
function code_toggle() {
if (code_show){
$('div.input').hide();
$('div.prompt').hide();
} else {
$('div.input').show();
$('div.prompt').show();
}
code_show = !code_show
}
$( document ).ready(code_toggle);
</script>
<form action="javascript:code_toggle()"><input type="submit" value="Code Toggle"></form>''')
from main import line_plot, levelset_plot, get_table_data, plot_order_of_accuracy
from itables import init_notebook_mode
init_notebook_mode(all_interactive=True)
line_plot(
data_name='free_energy_1',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Free Energy',
title="Free Energy, <i>r</i><sub>0</sub>=0.99<i>r</i><sup>*</sup>"
)
)
line_plot(
data_name='solid_fraction_1',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Solid Fraction',
title="Solid Fraction, <i>r</i><sub>0</sub>=0.99<i>r</i><sup>*</sup>"
)
)
line_plot(
data_name='free_energy_2',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Free Energy',
title="Free Energy, <i>r</i><sub>0</sub>=<i>r</i><sup>*</sup>"
)
)
line_plot(
data_name='solid_fraction_2',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Solid Fraction',
title="Solid Fraction, , <i>r</i><sub>0</sub>=<i>r</i><sup>*</sup>"
)
)
line_plot(
data_name='free_energy_3',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Free Energy',
title="Free Energy, <i>r</i><sub>0</sub>=1.01<i>r</i><sup>*</sup>"
)
)
line_plot(
data_name='solid_fraction_3',
benchmark_id='8a.1',
layout=dict(
x='time (s)',
y='Solid Fraction',
title="Solid Fraction, <i>r</i><sub>0</sub>=1.01<i>r</i><sup>*</sup>"
)
)
from main import levelset_plot
levelset_plot(
data_name='phase_field_1',
benchmark_id='8a.1',
layout=dict(
levelset=0.5,
range=[-10, 10],
title="Solid / Liquid Boundary at t=200, <i>r</i><sub>0</sub>=0.99<i>r</i><sup>*</sup>"
),
mask_func=lambda df: (0.1 < df.z) & (df.z < 0.9)
)
from main import levelset_plot
levelset_plot(
data_name='phase_field_2',
benchmark_id='8a.1',
layout=dict(
levelset=0.5,
range=[-10, 10],
title="Solid / Liquid Boundary at t=200, <i>r</i><sub>0</sub>=<i>r</i><sup>*</sup>"
),
mask_func=lambda df: (0.1 < df.z) & (df.z < 0.9)
)
from main import levelset_plot
levelset_plot(
data_name='phase_field_3',
benchmark_id='8a.1',
layout=dict(
levelset=0.5,
range=[-10, 10],
title="Solid / Liquid Boundary at t=200, <i>r</i><sub>0</sub>=1.01<i>r</i><sup>*</sup>"
),
mask_func=lambda df: (0.1 < df.z) & (df.z < 0.9)
)
plot_order_of_accuracy(
data_names=['phase_field_3', 'phase_field_4', 'phase_field_5'],
benchmark_id='8a.1',
keys=['x', 'y', 'z'],
rangex=[-50, 50],
rangey=[-50, 50],
nx=1000,
ny=1000,
layout=dict(
title='Order of Accuracy',
labels=dict(x=r'Δx', y=r'L<sup>2</sup> Norm')
)
)
Table of 8a.1 benchmark result uploads
get_table_data('8a.1')